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ABSTRACT 



Aims. We study the inner region (~1.0 AU up to a few 10 AUs) of the circumstellar disk around the "classical" TTauri star RYTau. 
Our aim is to find a physical description satisfying the available interferometric data, obtained with the mid-infrared interferometric 
instrument at the Very Large Telescope Interferometer, as well as the spectral energy distribution in the visible to millimeter wave- 
length range. We also compare the findings with the results of similar studies, including those of intermediate-mass Herbig Ae/Be 
stars. 

Methods. Our analysis is done within the framework of a passively heated circumstellar disk, which is optionally supplemented by 
the effects of accretion and an added envelope. To achieve a more consistent and realistic model, we used our continuum transfer code 
MC3D. In addition, we studied the shape of the 10 fim silicate emission feature in terms of the underlying dust population, both for 
single-dish and for interferometric measurements. 

Results. We show that a modestly flaring disk model with accretion can explain both the observed spectral energy distribution and 
the mid-infrared visibilities obtained with the mid-infrared infrared instrument. We found an interesting ambiguity: a circumstellar 
active disk model with an added envelope, and a lower accretion rate than in the active disk model without envelope, could represent 
the observations equally as well. This type of model with the envelope should be considered a viable alternative in future models of 
other TTauri stars. The approach of a disk with a puffed-up inner rim wall and the influence of a stellar companion is also discussed. 
We also investigate the influence of various fit parameters on the outcome of the radiative transfer modeling. From the study of the 
silicate emission feature we see evidence for dust evolution in a T Tauri star, with a decreasing fraction of small amorphous and an 
increasing fraction of crystalline particles closer to the star. 
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1. Introduction 

TTauri stars are known as precursors of low-mass main se- 
quence stars (^2-3 Mq). In contrast to main sequence stars, 
their characteristic properties are strong emission line radiation 
(e.g., Balmer a) and excessive continuum radiation observed in 
the UV, infrared and the millimeter (mm) wavelength range of 
their spectral energy distribution (SED). It has been shown that 
the spatial distribution of circumstellar dust in a disk or an enve- 
lope that is primarily exposed to stellar radiation, is responsible 
for the excess radiation in the infrared wavelength range (e.g., 
Adams et al. .l987) . while accretion of circumstellar material re- 
sults in the UV excess and strong emission line radiation (see 
Hartmannll998'for a review). 

The extraordinary interest in the inner region of a circumstel- 
lar disk results from the assumption that the formation of planets 
is favored there (see Nagasawa et al. 2 006, Wiinsch et al. 120051 
Klahr l2004l l. While mm observations probe cooler outer disk re- 
gions and layers close to the midplane of circumstellar disks, 
observations in the mid-infrared (MIR) wavelength regime are 
more sensitive to warmer (250 K ^ T ^ 1000 K; see Schegerer 
et al. 2006: Fig. 1) disk regions, such as the surface of the in- 
ner regions where dust is directly iiTadiated by the central star 
However, studies of the inner circumstellar regions (~1 AU) of 



objects in the closest star forming regions are only feasible by 
interferometric observations. 

Strong emission features in the MIR range at 10 fim and 
20 yum corresponding to the Si-O stretching and bending modes 
of silicate grains, are assumed to result from absorption and ree- 
mission processes in optically thin dust layers of the circumstel- 
lar disks. While silicate grains are expected to be initially amor- 
phous and small (^ 0.1 //m; Mathis et al. 119771 MRN there- 
afterfl the crystallization of amorphous silicates starts at tem- 
peratures of -1200 K (e.g., Gail |1998l l. Moreover, high dust 
densities and turbulent processes in the interior of circumstel- 
lar disks favor dust grain growth to dust pebbles (e.g., Blum et 
al. 2000; Johansen et al. 2005). The shape of the emitted silicate 
feature allows the estimation the predominant stage of the dust 
evolution in a young stellar object (YSO). Different degrees of 
crystallization and grain growth have already been shown in a 
large sample of T Tauri stars of different ages and stellar masses 
(e.g., Schegerer et al.|2006|l. As temperature and density increase 
in circumstellar disks with decreasing distance to the central star, 
crystallinity and grain size sensitively depend on the radial po- 
sition of the dust in a circumstellar system (e.g., Beckwith et 
al. 1999; WeidenschiUing 1997; Gail 2003 Fig. 28). In fact, ob- 
servations with MID-infrared Interferometric instrument (MIDI) 
have already revealed a correlation between the radial position 
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and the evolutionary stage of silicate dust in circumstellar disks 
around Herbig Ae/Be (HAeBe) stars, which are the more mas- 
sive counterparts of T Tauri stars (Leinert et al. 2004; van Boekel 
et al. l2004] i. 

The density, temperature and compositional structure of cir- 
cumstellar dust disks and surrounding envelopes have been 
the central issue of many former studies (e.g., Chiang & 
Goldreich [19971 D'Alessio et al. 2001). Different modehng ap- 
proaches have been tried to quantitatively explain and reproduce 
phenomenons like excess radiation, shapes of emission lines 
(e.g., Muzerolle et al. 2004; Natta et al. 12000 1. flux variations 
(e.g., Herbig et al. 2003 ) and intensity distributions of images 
(e.g., Lucas et al. 1997). However, the evolution of (inner) disk 
structure and its correlation to dust evolution is still unclear (e.g., 
Millan-Gabet et al. 2006 ; Beckwith et al. i 19991 ), and has been un- 
derestimated in actual modeling approaches, or mainly reserved 
for theoretical studies (e.g., Gail 1998 ). However, interferomet- 
ric observations in the MIR wavelength range, which are now 
available, are sensitive to the inner disk structure where warm 
dust dominates. Including their sensitivity for the silicate feature, 
the correlation between inner disk structure and grain evolution 
can be directly studied. 

In this paper we focus on modeling of the SED and spectrally 
resolved N band visibilities, which we obtained for the T Tauri 
star RY Tau with MIDI. The key questions of this study are the 
following: Is it possible to simultaneously model the SED and N 
band visibilities of RY Tau solely by an externally, i.e., passively 
heated, disk? Do different extensions of this model reproduce 
the observations, simultaneously? What do we learn about the 
(silicate) dust composition of the disk? 

The result of previous measurements of RYTau are pre- 
sented in Sect. |2l In Sect. [3] we outline the observations of 
RY Tau and the subsequent data reduction. We present the ra- 
diative transfer code and the basic dust set of our modeling ap- 
proach in Sect. H] In the following Sect. |5] we compare the re- 
sults of the different modeling approaches we used, i.e., the ative 
disk model with and without an envelope, and point to supple- 
ments. In Sect. |6] the dust composition of the upper disk lay- 
ers and its dependence on the radial distance from the central 
star is studied. Finally, in Sect. [T] we draw comparisons between 
the used models, refer to previous results, including the mod- 
els of Akeson et al. (120051 ) where near-infared (NIR) visibilites 
were modeled, and discuss discrepancies. Furthermore, we in- 
vestigate the possibility of the existence of a stellar companion 
and compare RYTau with HAeBe stars. Section [8] summarizes 
our results. 

2. Previous measurements 

RYTau, demonstrably observed for the first time in 1907 
(Pickering '1907^, is a well-known T Tauri star (Joy 119451 ) that 
belongs to the Taurus-Auriga molecular cloud at a distance of 
~ 140 pc. Table[T]shows the main properties of this object, which 
are obtained from previous measurements. Photometric fluxes 
are listed in Table |2l 

RYTau is a UXOri-type star, i.e., this T Tauri object has re- 
vealed irregular photometric variability in the visible and NIR 
wavelength range. During several months in 1983/84 and in 
1996/97, its visible brightness increased from ~llth to ~9th 
magnitude and decreased again to its initial value (Herbst & 
Stine [19841 Zajtseva et al. 1985"; Petrov et al. 1999'; Herbst & 
Shevchenko[i999 ). Such a rare but strong variability is conven- 
tionally explained by variations of the obscuration of the central 
star caused by an inclined circumstellar disk and an envelope 



Table 1. Observed properties of RYTau. References - 1: 
Perryman et al. (1 9971 ); 2: Bertout et al. (T^9); 3: Calve t et 
al. (2004); 4: Mora et al. (|2001); 5 : Beck with et al. ([T990l ); 6: 
Akeson et al. (120051 ); 7: Siess et al. ([T9991 ) 



Parameter Value Reference 



RA (J2000.0) 


04 21 57.4 


1 


DEC (J2000.0) 


+28 26 36 


1 


Distance 


134!^4pc 


2 


Visual Extinction 


(2.2 ± 0.2) mag 


3 


Spectral type 


F8 III 


4 


Stellar mass 


1.69 Mo 


5 


Stellar luminosity 


12.8 Lo 


6 


Accretion rate 


7.8 X 10-** Mo/yr 


3 


Age 


(6.5 ± 0.9) Myr 


7 



Table 2. Photometric measurements of RY Tau. The data of dif- 
ferent measurements are averaged and the standard deviation is 
determined. For convience, all fluxes are given in Jansky and 
in magnitude. Conversion factors are taken from Leinert (19971 ). 
References - 1: Calvet et al. (119971 ); 2: Gezari-Catalog (Gezari 
et al. il999D; 3: 2MASS-Catalog (Cutri et al. 2003 ); 4: Rydgren 
et al. ([T9761 ); 5: EHas (1978 ); 6: Ha nner et al. (1998j; 7: Strom 
et al. ([T9891 ); 8: Harvey et al. ( [19791 ); 9: Mora et al. ( [2U0T1 ) 



wavelength [/im] flux [Jy] ([mag]) Reference 

0.36 (U) 0.04 + 0.01 (11.72 + 0.16) 1 

0.45(B) 0.14 + 0.01 (11.23 + 0.08) 1 

0.55 (V) 0.33 + 0.01 (10.12 + 0.04) 1 

0.64 (Rc) 0.59 + 0.03 (9.3 + 0.05) 1 

0.79 (Ic) 1 .0 + 0.06 (8.52 + 0.06) 1 



1.25 (J) 1.6 + 0.8(7.64 + 0.40) 2,3 

I . 65 (H) 2.3 + 0.9 (6.48 + 0.30) 2, 3 
2.20 (K) 3.8 + 0.5(5.56 + 0.13) 2,3 
4.80 (M) 6.6 + 2.0 (3.4 + 0.28) 4, 5 

II. O(N) 20.0 + 0.3 6 
25 28 + 3 7 
60 18 + 4 7 
100 12 + 5 7,8 
1300 0.23 + 0.02 9 



(e.g. Eiroa et al. [20021 ). Smaller variations (AV ^ 0.1, A7 » 0.2, 
AK ss 0.2) in the range of several days were also detected (Eiroa 
et al. 120021 ). By a comparison between the maximum and min- 
imum brightness of the object the photometric measurements, 
listed in Table |2| and our observations with MIDI coiTespond to 
the "quiescent" state of the object, i.e., close to the photometric 
minimum. 

There is a wide range of values measured for the vi- 
sual extinction Ay of RYTau (Kuhi 1974: 1.3 mag; Cohen & 
Kuhi il979i (1.9+0.2) mag; Strom et al. il989i 0.6 m ag; Beckwith 
et aL [19901 2.7 mag; Kenyon & Hartmann [T9951 1.8 m ag). We 
adopt a value recently derived by Calvet et al. (119971 ): Ay = 
(2.2 ± 0.2) mag. The level of veiling in the visible range of the 
spectrum is low (^ 0.1; Basri et al. |1991| Hartigan et al. 119951 
Petrov et al. '1999^ but markedly higher in the infrared range 
(> 0, 8 ; Folha & Emerson 1999 ). 

A potential duplicity/multiplicity was not found by Leinert et 
al. (11993 ) by using NIR speckle interferometry reaching a spa- 
tial resolution between 0.13" and 13" but the regular variation 
of the photocenter, found by the astrometric measurements of 
HIPPARCOS, could be a hint for a companion with a projected 
minimum distance of 3.27 AU (23.6 mas) and a position angle 
of 304° + 34° (Bertout et al. [T9^ . 
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3. Interferometric observations and data reduction 

3.1. Observing sequence 

RY Tau was observed with MIDI/VLTI (Very Large Telescope 
Interferometer; Leinert et al. '2003a) in 2004, November P' and 
4*, within the scope of guaranteed time observations. The dates 
and universal times (UT) of the observations, as well as sky- 
projected baseline lengths (L) and position angles (PA) of the 
interferometer, are listed in Table[3] An observing sequence con- 
sists of the following steps: 



Single-telescope imaging is used for a highly precise acqui- 
sition of the object within a field-of-view of 2 x 2" in order 
to guarantee a maximum overlap of the telescope beams. 
After the beam combiner is introduced into the optical 
path, the characteristic interference pattern (fringe pattern) 
is found around an optical path difference (OPD) of zero. A 
low resolution prism (A/ 6 A » 30), which is put in the optical 
path of the combined beams, allows us to obtain spectrally 
resolved interferograms and the wavelength-dependent cor- 
related flux F^on-i^). 

The spectrum Fiof^i(A) is determined by single-telescope ex- 
posures that are recorded on the same detector pixels as the 
fringe signal Fcon(/l)- By definition, the spectrally resolved 
visibility V(A) is obtained as the ratio of the correlated and 
the total flux 



111. 



V{A)^ 



^'total('i) 



(1) 



This visibiUty is still biased by instrumental eff'ects. 
The transfer function of the instrument is determined by 
the observation of calibrator stars before and after the ob- 
servation of the scientific target. A known transfer function 
is required for the elimination of the instrumental and at- 
mospheric influences. The calibrators are selected for their 
known diameter, the absence of strong photometric variabil- 
ities and companions, and their angular distance to the sci- 
entific target. Some of the calibrators were also used for ab- 
solute flux calibration. The error of our calibrated visibility 
V(A) is the 1 cr deviation that is obtained by the observations 
of dififerent calibrators in one night (see TableO. 



Table 3. Summary of the MIDI observations of RY Tau and cal- 
ibrators. The dates, UT, L (in m) and PA (in degrees, measured 
from North to East) of the sky-projected baselines are listed. The 
airmass AM in the right column is given for the time of fringe 
tracking. The observations with a projected baseline of 79 m and 
81m provided an almost identical result (see Fig.[T]i. 



Instrumental informations and observing procedures are 
described by Leinert et al. ( |2003bl |2003cl |2004b and 
Ratzka (2005). 

3.2. Data reduction 

The reduction procedure of MIDI data is complex and has been 
described in detail by Leinert et al. (2 004) , Ratzka (120051 ). and 
Jaffe (120041) . The data obtained with MIDI were reduced with the 
MIA software that is based on power spectrum analysis and the 
results were independently confirmed by using the EWS soft- 
ware. The EWS software contains a coherent integration algo- 
rithm, which involves a kind of shift-and-add in the complex 
plane. Both reduction software packages are publicly availableQ 

3.3. Observationai resuits 

The resulting wavelength dependent visibility curves for the 
three baselines are shown in Fig. [1] including 1 cr error bars. 
The spectrophotometry of the silicate emission band /^totai('^), 
obtained during the measurements, is shown in Fig.]?] and also 
included in the SED of RY Tau (Fig. |6]l. Figure [T] also shows 
the observed correlated fluxes Fcoit(^), formally obtained as a 
product of these two quantities (Eq. [T]). We refer to the near- 
coincidence of the observations at baseline lengths of 79 m and 
81 m. 




* bun] 



Fig. 1. The spectrally resolved, calibrated visibility data derived 
from our MIDI observations (Table|3]l. The error bars are the 1 cr 
deviations that result from the observations of diff'erent calibra- 
tors in one night. 



Date 


UT 


Object 


L 


PA 


AM 


Nov. V\ 


2004 


3:54 


-4:16 


HD 25604 


74 


96 


1,68 


Nov. P', 


2004 


4:37 


-4:56 


RYTau 


79 


97 


1.83 


Nov. P', 


2004 


4:58 


-5:07 


RYTau 


81 


95 


1.79 


Nov. P', 


2004 


5:54 


-6:10 


HD49161 


64 


87 


1.56 


Nov. P', 


2004 


7:03 


-7:17 


HD 31421 


89 


82 


1.28 


Nov. 4'\ 


2004 


0:01 


-0:26 


HD 178345 


57 


146 


1.42 


Nov. 4*, 


2004 


2:19 


- 2:47 


HD 188603 


46 


169 


2.47 


Nov. 4'\ 


2004 


3:16 


-3:54 


HD 25604 


61 


117 


1.75 


Nov. 4'\ 


2004 


5:11 


-5:29 


HD 20644 


59 


102 


1.69 


Nov. 4*, 


2004 


7:18 


-7:36 


HD 37160 


61 


107 


1.21 


Nov. 4'\ 


2004 


7:44 


- 8:00 


RYTau 


49 


92 


1.95 


Nov. 4'\ 


2004 


9:00 


-9:23 


HD 50778 


61 


113 


1.04 



4. Tools 

4.1. MC3D - Monte Carlo code for radiative transfer 

In contrast to many previous investigations where the radiative 
transfer function has been solved (e. g., Sonnhalter et al. 119951 
Chiang & Goldi-eich |1997i Dullemond et al. 120011) . we use the 
well-tested code MC3D, which is based on the Monte-Carlo 
method (Wolf et al l 19991 Pascucci et al. 120041 ). Considering an 



^ |http : //www ■ mpi a-hd ■ mpg . de/MIDISOFT/ and 

I http: //www. strw. leidenuniv.nl/~koehler/MIA+EWS-Manual/| 
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axially symmetric object, we assume a two-dimensional geome- 
try in a polar coordinate system {r,0) with a logarithmic grid spa- 
cing in r and a uniform grid spacing in 9. Heating sources like the 
central star, accretion effects, and heated dust grains determine 
the temperature distribution. The flux of the central star is deter- 
mined by the theoretical stellar atmosphere model provided by 
the Kurucz (1994). The product of the dust-specific absorption 
efficiency 2ab,s(^), the grain surface and the blackbody emission 
Bv(Tdust) represent the flux that a dust grain with a temperature 
of Tdust reemits in local thermal equilibrium (Kirchhoff's law 
of local thermal equilibrium). Sources like the central star and 
accretion are treated as blackbody emitters. Gas molecules and 
atoms are not considered in our models. 

The radial density distribution of the disk is given by the 
surface density profile 

2(r) = Eo ■ (2) 

with the radial coordinate r, an exponent p and a constant 
So = 100 gcm~^ ( Weidenschilling ' 1 977i The vertical density 
distribution is calculated self-consistently assuming hydrostatic 
equilibrium, i.e., the balance of gravitational and thermal pres- 
sure (Schegerer et al., in prep.). 

After temperature and density distribution have been deter- 
mined, the SED and the projected image of the star and its cir- 
cumstellar environment, considering an inclination angle j?, are 
calculated. The resolution of the image is by a factor of ~10 bet- 
ter than the resolution of the observations. For projected base- 
lines of 79 m and 49 m, our observations with MIDI reached 
spatial resolutions of ~1.8 AU and ~2.8 AU at a distance of 
~14Opc0 

4.2. Dust model 

The infrared excess that is emitted from YSOs originates from 
heated dust in the circumstellar environment. Assuming com- 
pact, homogeneous, and spherical dust grains, their optical prop- 
erties, such as scattering and extinction cross sections, are deter- 
mined by Mie scattering-theory from the measured complex re- 
fractive index of the specific material (Bohren & HufFman ll983l l. 
In our modeling approach, we assume a dust mixture of "astro- 
nomical silicate" and graphite with the relative abundances of 
62.5% for astronomical silicate and 37.5% for graphite (Draine 
& Malhotra 1993). The dielectric function of astronomical sili- 
cate was formerly synthesized by Draine & Lee (1 19841 1 in order 
to reproduce the extinction of different silicate compounds in 
interstellar space. We consider an improved version of this di- 
electric function (Weingartner & Draine 20000, which was re- 
cently confirmed by a study of the interstellar extinction in the 
NIR wavelength range (Indebetouw et al. 120051 ). For graphite we 
adopt the | : | ratio with Ky - [Ky{e\\) + 2Ky(e±)]/3, where e\\ and 
e± are the components of the graphite dielectric tensor for the 
electric field parallel and perpendicular to the crystallographic 
c-axis and Ky is the mass absorption coefficient. The strongly ab- 
sorbing graphite grains efficiently contribute to the heating of the 
dusty circumstellar environment. The ratio of the extinction ef- 
ficiency factor of carbon dust to the extinction efficiency factor 
of silicate dust is -10 in the NIR wavelength range (Draine & 
Lee [T984l Wolf & Hillenbrand|2003). 



' The spatial resolution power R of an interferometer is given by the 
ratio of the observing wavelength A and the sky-projected distance be- 
tween a telescope pair, i.e., the effective baseline length B. R = A/(2B). 
See ,http : //www . astro . princeton . ed u/ -draine. 



We consider a grain size power law n{d) oc oT^-'' with 
fl„in ^ fl ^ flmax> whcrc «(fl) is the number of dust particles with 
radius a. This power law was formerly found by MRN studying 
extinction of interstellar carbon and silicate with typical sizes be- 
tween 0.005yum - O.Oljum and 0.025 ^im - 0.25jum, respectively. 
This grain size power law has akeady been used in former mod- 
eling aproaches of YSOs. We use a minimum particle size of 
flmin = 0.005 /im in all of our modeling approaches. 

The maximum grain size cimax affects the mass absorption 
coefficient Ky{ci) of dust, i.e., Ky{d) increases with a^ax for sizes 
up to a few mm. The mm flux depends on the disk mass Mdisk 
and the absorption coefficient Ky(d), i.e., Omax, in particular, when 
assuming an optically thin disk in the mm wavelength range. 
Correspondingly, the spectral index a (Fy oc v") decreases 
with an increase of cjmax from an absolute value of ~4 (for 
flmax < 0.1 ;um) to ~2 (only for large bodies). The correlation be- 
tween the spectral index a and the maximum dust grain size Omax 
was formerly studied by D' Alessio et al. (.2001) , while Wood et 
al. (120021 ) also investigated the correlation between disk mass 
and mm flux. 

In our modeling approach, we have found that a dust distri- 
bution with a maximum grain size of Omax = 0.25/im and the 
above mentioned grains size power law for silicate and carbon 
generally underestimates the mm flux unless a disk mass in the 
range of ~ 1 .0 M© is assumed. Circumstellar disks with such high 
masses are potentially gravitationally unstable (e.g., Laughlin 
& Bodenheimer [1991 Boss llOOOl Lodato & Berlin |250TT ). 
Moreover, an upper grain size of flmax - 0.25 jum results in 
a too steep mm slope, in contrast to the measured spectral in- 
dex a. Using the Very Large Array for their mm measurements, 
Rodmann et al. (120061) found a spectral index of a = 2.55 + 0.09 
for RY Tau and derived a maximum grain size of Omax = 1 mm. 
The latter maximum grain size is used in our modeling ap- 
proach. We have to mention that the spectral index a provides 
only a lower limit for the maximum dust size as it converges for 
flmax > 1 mm. 

Although the mm wavelength range of the observed SED can 
sufficiently be simulated considering grain sizes up to 1 mm and 
relatively low disk masses (< 10"^ Mq), the spectral contribution 
in the NIR wavelength range strongly decreases with increasing 
maximum grain size. The dust particles with Omax > 1 /^m can be 
less effectively heated than the smaller particles of the canonical 
MRN distributionQ This effect is the reason why we implement 
a two-layer dust model in our modeling approach. The disk in- 
terior contains a maximum dust grain size of 1 mm while the 
MRN grain size distribution with fl^ax - 0.025 pccs\ is used in 
the upper disk layers where the optical depth tn in N band, mea- 
sured vertical to the disk midplane, falls below unity. Such a 
division of the disk is based on the idea of the favored settling 
of larger dust grains. Furthermore, dust particles are assumed to 
mainly grow in the denser regions of the disk close to the mid- 
plane (e.g., Lissauer |1993| Blum & Wurm |2000] ). Similar disk 
models with two or more different dust layers have already been 
proposed by Chiang & Goldreich ( 1997 ) and used by Whitney et 
al. ( 120031 ), for instance. However, it is still an open question how 
strongly dust grains are mixed in the circumstellar environment 
(e.g., D' Ale ssio et al. [T997l Gail [T998l McCabe et al. 2003; Wolf 
etal. 120031 ). 

In order not to determine the temperature distribution of 
each single dust component and to accelerate the radiative trans- 
fer simulations, we construct an "artificial" particle with optical 



' For all different dust sets we assume a constant exponent of -3.5 
for the grain size power law n{d). 
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constants that are derived by averaging the optical properties of 
carbon and astronomical silicate of different sizes in each dust 
layer. Such an approach was justified by Wolf (.2003 j . 

5. Models of the density structure 

In the following subsections we will present our approaches to 
model the SED and the MIR visibilities that we obtained with 
MIDI. An active disk model is our favorite approach. Such a 
model was also used to reproduce SED and K Band visibili- 
ties of RY Tau obtained with the Palomar Testbed Interferometer 
(Akeson et al. 120051 ). Additionally, an active disk model with an 
envelope and potential supplements are discussed. One of the 
main issues of this study is to clarify if the models we use can be 
distinguished. 

The existence of a passively-heated circumstellar disk in our 
analysis is beyond all questions. The paradigm that the forma- 
tion of a disk is one of the evolutionary stages of circumstel- 
lar stuctures has been finally confirmed by the images of many 
YSOs in different wavelength ranges (e.g., Padgett et al. 119991 
Mannings & S argent 79971 l2000b . The TTauri star RYTau is a 
Class II object (Kenyon & Hartmann[l995) where a surrounding 
disk has already formed. Therefore, a passively heated disk is a 
basic ingredient and will be retained in our different modeling 
approaches. 

The disk model is characterized by the disk mass Mdisifl 
the inner Ri^ and outer disk radius Ront, the inclination angle &, 
and an exponent p, which represents the radial dependence of 
the surface density S (see Eq. |2]i. The inner radius Ri^ is given 
in advance, but can also be considered as a starting value that 
is iterated until the temperature at the inner radius falls below 
the the sublimation temperature of 1600 K (Duschl et al. 1 1 9961 ). 
Properties of the central star, like the stellar temperature T*, stel- 
lar luminosity L* and stellar mass M*, are additional model pa- 
rameters but these quantities are well constrained by previous 
studies (see Table[T]i. 

5. 1 . Active disk 

Our active disk model is a passively heated (dust) disk where ac- 
cretion effects are added. The existing MC3D radiative transfer 
code therefore has to be extended. The implemented accretion 
effects are briefly described in Appendix lAl 

Apart from the parameters of the disk and the star, our ac- 
cretion model requires three additional model parameters: ac- 
cretion rate M, boundary temperature Tbnd of the accreting re- 
gions on the surface of the star, and a magnetic truncation radius 
Rhnd- Defining the inner radius of the gaseous disk inside Ri„, 
Konigl ( 119911 showed that the truncation radius is not an inde- 
pendent quantity, but depends on stellar radius, mass, accretion 
rate and magnetic field strength. However, as we do not know 
the exact magnetic field strength of RY Tau, we fix the boundary 
temperature and truncation radius to 8000 K and 5 Ri,, respec- 
tively. Both values were already used for the same object by the 
study of Akeson et al. (2005') and were justified by the assumed 
large magnetic field of RY Tau in the range of a few kilo-Gauss. 

The best fit parameters for this model are given in Table |4] 
and the model is compared to the observations in Fig. |6l The 
accretion luminosity is 1.2 Lq. This model suffers from the fol- 
lowing deficiency: the far-infrared (FIR) wavelength range in 
the SED is slightly overestimated. A potential improvement of 



* We assume a gas-to-dust mass ratio of 100: 1 . 



this active disk model could be a "truncated outer disk", for- 
merly suggested by Lucas & Roche (1997 ) and recently used by 
Preibisch et al. ( 120061 ). For this, the primary density distribution 
of the disk is truncated at an outer radius /?td by multiplying the 
surface density E(r) (Eq. |2]) with a Fermi-type function. With a 
constant Ctd > 

^TD(r) = -— — — — (3) 

1 -Hexp(CTD^ - Ctd) 

This decreases the extended IR emission while the inner disk 
structure is modified much less than the outer disk regions. It 
could have its physical origin in a stellar companion which trun- 
cates the outer disk region. In fact, Bertout et al. ( 119991 ) pointed 
to a potential binarity of the system. Our study has shown that 
only a truncation radius of the order of 10 AU could effectively 
decrease the computed SED (Preibisch et al. 2006). However, 
this truncation radius is still too large to correspond to the find- 
ings of Bertout et al. (119991 ). Moreover, such a truncation radius 
is so small that the MIR visibilities increase. Another possibil- 
ity to decrease the extended IR emission is dust settling that re- 
sults in a flattening of the flared disk as suggested by Miyake & 
Nakagawa (1995) and Dullemond & Dominik (2006). As outer 
disk regions are less heated, a less effective flared disk would 
also decrease the visibility. 

Studying SED and visibility, the inclination ■& and the posi- 
tion angle PA of the object cannot be clearly derived. Two visi- 
bility points are not sufficient to derive these values. With respect 
to the SED only an upper limit of ^65° can be determined. This 
corresponds to the angle of the line of sight where the optical 
depth Ty in the model exceeds unity. 

The outer disk radius 7?out - 270 AU of this model is 
larger than the results from previous measurements. Andrews & 
Williams (2006) derived an outer disk radius of 150 AU with the 
Submillimeter Array (SMA), while Rodmann et al. (2006 ) found 
an outer disk radius of only (90 ±15) AU with 7 -mm continuum 
observations at the Very Large Array (VLA). Fig. |2] shows at 
which stellar distances the MIR flux arises in our model. This 
radial flux distribution confirms that it has its origin in the in- 
ner disk regions (< 10 AU), mainly, that can only be observed 
with interferometric methods. But it also shows that MIDI is not 
sensitive to outer disk regions (> 40 AU), including the outer 
disk radius. A modeling of mm maps of this object allows to 
determine the outer disk radius but this is out of the scope of 
this paper Additionally, it should be underlined that the distri- 
bution of the NIR flux from the disk culminates at ~1 AU for 
this model. 

The disk mass of this model Mdisk = 4 x lO"-' M© is by 
a factor of 3 smaller than the value that was found by Akeson 
et al. ( 120051 ). In this context we have to underline that the disk 
mass strongly depends on the dust set that is used in the model 
(see Sect. 14.2b . The Fig. [3] left, shows the SED for our active 
disk model but with the MRN dust set (Omax = 0.25 yum), only. 
All other model parameters are adapted. In this figure it can be 
seen that the resulting flux in the mm wavelength range strongly 
declines for a^ax = 0.25 jim. 

Finally, we attempted to reproduce SED and N band visibil- 
ities with a passive disk model, solely, instead of implementing 
accretion effects, additionally. Such passive disk models for dif- 
ferent parameter sets were calculated but without reproducing 
the SED and MIR visibilities, sufficiently. Each modeled SED 
of these simulations suffers from similar deficiencies (see upper 
panel in Fig.|6]): in contrast to the photometric measurements the 
predicted NIR flux between A ^ 3 fim and 8 yum is generally 
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Fig. 3. Left: Deficiency of the active disk model at long wavelength if only dust with standard size MRN distribution is assumed 
(<5imax = 0.25 /im; compare middle panel in Fig.|6]l. Observations are represented by crosses with error bars (s. Table|2]i. The predicted 
slope into the mm range is too steep. The dashed line shows the unreddened stellar flux. Right: The same as the left figure but an 
maximum dust grain radius of Omax = 0.25 fj.m is assumed. 




rodius [AU] 

Fig. 2. Radial flux distribution A//Ar for a wavelength of 8 fj.m 
(solid), 10 yum (dotted), 12 yum (dashed curve) and 2.4 /im 
(dashed-dotted curve). The mean flux /, emitted at a radius r, 
is multiplied with Anr^ and normalized by the maximum of the 
radial flux distribution / at a wavelength of A - E fim. 



underestimated. Moreover, the model is too strongly spatially 
resolved (visibilities which are too low) in comparison with the 
measured visibilities. 

This lack of NIR radiation in "naked" passive disk models 
was reported by Hartmann et al. in 119931 Obviously, accretion 
can generate the missing NIR excess in a region which cannot 
resolved by the interferometer. In the following subsections we 
present two further modifications (the puffed-up inner rim wall, 
an envelope) which have been considered to reproduce distinct 
NIR excess. 

5.2. Inner rim wall 

In addition to an accretion model that is described in Sect. 15.11 
Natta et al. (2001 ) suggested that a puffed-up wall at the inner 
disk rim could also produce an increased NIR excess. Dust at the 
inner edge of the disk is strongly heated by direct stellar irradia- 



tion. Based on the idea of hydrostatical equilibrium this heating 
of the inner disk edge is assumed to cause an expansion of the 
dust layers perpendicular to the midplane of the disk. DuUemond 
et al. ( 120011) established an analytical model of the puffed-up, 
inner rim wall as a supplement to the Chiang-Goldreich model 
(Chiang & Goldreich |1997l l and defined the "inner rim scale- 
height" i/rim as follows: 

^^rim = A'rim/iin, (4) 

where ;^'iin, is a dimensionless factor, larger than unity. The quan- 
tity /zin represents the scale-height of the inner disk without 
extra vertical expansion^ Simulations of young stellar objects 
(YSOs) where the inner rim wall has been successfully applied, 
have been published by Dominik et al. (|2003 ), Pontoppidan et 
al. (.2005 ) . and Eisner et al. (2005ill, for instance. The reproduc- 
tion of interferometric observations in the NIR regime assuming 
simple, analytical ring models (Millan-Gabet et al. '200 H has 
also been used as a further confirmation of the inner rim wall 
(Natta et al. 2001). Furthermore, a puffed-up inner rim, shad- 
ows adjacent areas of the disk from direct stellar radiation ("self- 
shadowed disk"). This effect is used to explain the spectral shape 
of the FIR excess from YSOs and to explain the classification of 
Group I and Group II sources (Meeus et al. 2001 ; Dullemond & 
Dominik 2004a). Arising from a flared disk a Group I source has 
a flat SED over the entire infrared wavelength range. The SED of 
a Group II source, however, strongly declines towards the FIR. 
Such a decline was explained by a self-shadowed disk. 

The puffed-up inner rim wall is still a controversial topic 
(e.g., Millan-Gabet et al. 12006). In particular, it was shown 
that a puffed-up inner rim wall does not generally emit enough 
radiation to cause the observed NIR excess in contrast to an 
envelope similar to the one that we implement in Sect. 15.31 
(Vinkovic et al. I2006I I. Another open issue is the static sta- 
bility of the proposed sharp, inner rim wall (Dullemond et 
al. i2001l l. Considering the sublimation temperature of the used 
dust species, Isella & Natta (2005 ) revised the previous inner rim 
model by a more rounded-off inner rim. Monnier et al. 2006 ob- 
served RY Tau in the NIR wavelength range using the Infrared 

' The scale-height is defined as the vertical distance from the mid- 
plane where the density has decreased by a factor e x 2.718 (Euler's 
constant). 
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Optical Telescope Array (IOTA) where a spatial resolution of 
^1 AU has been reached, comparable to our MIDI observations. 
In fact, their modeling results were incompatible with the mod- 
els possessing vertical inner walls. 

In our model the vertical density distribution is calculated 
assuming hydrostatic equilibrium. Therefore, the potential for- 
mation of a puffed-up inner wall is included in a natural way. 
However, it requires special computational care. In order to de- 
tect the effect of the puffed-up inner rim wall the size of grid 
cells in the inner region of the disk model should be small 
enough. A too coarsely meshed grid results in too low, aver- 
aged cell temperatures and in the absence of a potential puffed- 
up inner rim wall. We use a polar coordinate system (r, ff) in 
our two-dimensional model with uniform steps in 6 (AO - 1°) 
and a logarithmic scale for r (most inner step Ar w 0.01 AU). 
Therefore, the inner grid cells have an approximate size of 
~(0.01 X 0.005) AU^. According to Dullemond et al. (l200TT l 
a puffed-up scaleheight of 0.05 AU up to 0.10 AU can be ex- 
pected. Figure|4]shows the scale-height of our active disk model 
of Sect. O 

The fact that we do not see an excessively puffed-up inner 
wall in our approach does not necessarily exclude such a phe- 
nomenon. It is possible that the T Tauri star RY Tau is still too 
faint. In their study Dullemond & Dominik (2004a) assumed 
HAeBe stars with temperatures in the range of Tetj = 10, 000 K 
and luminosities up to several 10 solar luminosities (see also 
Dominik et al. I2003I I. The optical depth of the inner disk re- 
gion also affects the formation of a rim wall. If the optical depth 
is small at the inner disk region, the radiation is not absorbed 
mainly at the inner rim of the disk but on a larger scale. In this 
context Dullemond & Dominik ('2004a) showed that the verti- 
cal height of the inner rim wall is notably boosted for exponents 
p ^ 4 for the surface density (Eq. |2]i. Finally, the properties of 
the dust grains that are used in the modeling approach can also 
affect the inner rim as Vinkovic et al. (2006) mentioned. In Fig.|2] 
it is shown that the inner rim of our model emits only a smaller 
fraction of the NIR radiation. 

Our simulations are based on a Monte Carlo approach. Its 
drawback is the presence of statistical noise in the result which 
could in principle blur an effect such as the puffed-up inner 
rim wall. In order to quantify an upper limit of the scaleheight 
of a potential puffed-up inner rim, we assume that the scale- 
height at the innermost 0.3 - 3 AU can be described by a 
quadratic polynomial {h oc r^). This polynomial is then fitted 
to the scaleheight. The corresponding standard deviation itsd 
between the quadratic function and the scaleheight represents 
the maximum height of a potential puffed-up inner rim. We get 
CSD ~ 0.002 AU which is much smaller than the value derived 
by Dullemond et al. ( l200TT l. 

The strength of the inner rim wall is still an open issue. It de- 
pends strongly on dust properties, radial disk structure and stellar 
properties that are used in the models. Future, highly spatially 
resolved observations in the NIR wavelength regime, that are 
sensitive to the hot inner edge of the disk, will us allow to decide 
to what extent the effect of a puffed-up inner disk wall exists and 
which observational effects are provoked by this phenomenon. 

5.3. Active disk witli a dusty envelope 

Many studies (e.g., Hartmann et al. '1993' Calvet et al. 1 19971 1 
have shown that a dusty envelope around YSOs and Class-I 
sources, in particular, substantially contributes to the observed 
NIR excess. With respect to its optical depth a dusty envelope 
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Fig. 4. Scaleheight versus disk radius of our disk model from 
Sect. 15.11 We do not see a puffed-up inner rim in the sense of a 
local maximum of the scaleheight at the inner rim and a follow- 
ing local minimum at slightly larger radii although the rim of our 
model catches a large fraction of stellar radiation. 

could even dim the stellar radiation in the visible wavelength 
regime. 

There are also several studies which justify an envelope 
structure around the star and disk of RY Tau. Vink et al. (.20031) 
found that changes of the polarization across the Ha Une of 
RY Tau are based on scattering effects due to an extended dusty 
envelope. Direct evidence for such a circumstellar halo has 
been provided by R and I coronagraphic, large-scale images of 
RY Tau (Nakajima & Golimowski 1995 ) , in particular, and NIR, 
scattered Ught images around different YSOs (e.g., Padgett et 
al. 1 19991 Allen et al. |2002| l, in general. Certainly, one challenge 
for interferometric studies is to decide whether the observed 
large-scale halo around RYTau extends down to the inner disk 
region. 

In the context of axisymmetric accretion models 
UMch (1976) created an infall model of circumstellar gas 
and dust in an envelope structure in order to reproduce the 
emission-line Ha and H/3 profiles of type I/II P Cygni objects. 
This ansatz has been successfully used in modeling infrared 
images of Class-I objects (e.g., Lucas & Roche 119971 Wolf et 
al. 120031 ). In contrast to UMch's approach we add a more simple 
spherical dust configuration to the disk model. The spherical 
envelope in our model is geometrically constrained by the inner 
(Rin) and outer (Rom) disk radius. With the density distributions 
of the disk pdisk('", 6), of the envelope Penv(f, 0), and the the 
position vector r as well as the coordinates r and we define 

Penwir, 0) = Ci ■ Pdisk(^in, 0) • j (5) 

where ci << 1 and C2 > 0. The constraint ci << 1 guaran- 
tees a low optical depth of the envelope and the possibility to 
observe the innermost region of the disk. Disk and envelope are 
combined by 

p(r) = Pdisk('-) for PenAr) ^ Pdiskir) 

and (6) 
P(r) = Pem(r) for Pe„y(r) > pdi.kir). 

This assumption ensures a smooth transition from disk to enve- 
lope. In this simple envelope -i- active disk model we do not im- 
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plement bipolar cavities although there are hints that such cavi- 
ties, which are caused by collimated outflows, i.e. jets, generally 
exist in YSOs (Edwards et al. 1993 ). Actually, the number of free 
modeling parameters should be as small as possible to ensure 
that we do not overdetermine our approach. Furthermore, based 
on images of YSOs at a similar evolutionary state as RY Tau, 
Eisner et al. (|2005 b^ mentioned that cavities have effects on the 
structure of scattered light emission, which has its maximum in 
the NIR wavelength range. 

The Figure |6] and Table |4] show the result and parameter set 
of our best envelope + active disk models. The accretion rate of 
the model M = 2.5 x 10"^ M^yr^ is by a factor of ~4 smaller 
than the accretion rate assumed in the previous model without 
envelope. However, both models reproduce the SED and the vis- 
ibilities. This result shows that both accretion and an envelope, 
have the same effects on the SED and visibilities in the NIR and 
MIR wavelength range. Moreover, we have to mention that the 
measurements could be reproduced without considering any ac- 
cretion effects. A comparison of these results follows in Sect.]?] 
The accretion luminosity is 0.3 Lq. In the visual range and for 
inclinations § ^ 65° the model is optically thin as can be seen in 
Fig. |5] However the envelope evokes an observational effect on 
the SED and MIR-visibilities. 




50 60 ^ 70 

inclination "d [] 

Fig. 5. Optical depth for a inclination of § in the optical 
wavelength range. For smaller inclinations the optical depth is 
smaller than unity ensuring the observations of the inner regions. 
However, this envelope still has effects on the SED and MIR vis- 
ibilities. 



6. Radial gradients of the dust composition in 
circumstellar disks 



As mentioned in Schegerer et al. (1200 61) a more advanced disk 
model considers the specific dust composition of the correspond- 
ing object instead of using a canonical MRN dust set with av- 
eraged optical quantities (see Sect. I4.2l i. In such dust models 
the following dust components are generally taken into account: 
carbon, which mainly contributes to the underlying continuum 
as well as amorphous and crystalline silicate dust which gen- 
erates the silicate features at ~10 fim and ~20 jim. Other dust 
species such as water ice or, more precisely, water ice mantles 
around carbon/silicate grains can also be considered (see Fig. 1 
in Chiang et al. |200TT l. 



First comparisons between the silicate feature measured in 
laboratory experiments and observationally-based silicate spec- 
tra of YSOs have been drawn by Jager et al. ( 1994). A widely 
accepted analysis to determine the silicate composition of cir- 
cumstellar dust is a -fitting method that was established by 
Bouwman et al. (120011 ). They assumed that the silicate emission 
feature has its origin in the optically thin surface layer of the 
circumstellar disk where it results from a linear combination of 
mass absorption coefficient (emissivity) a:, of different dust com- 
ponents /: 



F(v) = B(v, T) 



Co H- ^ C//c,(v) 

1=1 



(7) 



where Co and C, are fitting parameters which reflect the mass 
contribution of each component /. The quantity F{v) is the spec- 
tral flux at frequency v, /c,(v) represents the frequency-dependent 
mass absorption coefficient for dust component / and By{T) is 
the Planck function corresponding to a blackbody temperature 
T. As basic dust set for their -fitting routine for TTauri ob- 
jects Schegerer et al. ( 12006 1) used the following silicate species: 
small (0.1 jum) and large (1.5 jjm) grains of amorphous olivine, 
and pyroxene, as well as crystalline species such as forsterite, 
enstatite, and quartz. 

The single-dish spectrum Fy is measured with a single tele- 
scope. For each interferometric observation of RY Tau we ob- 
tained a correlated spectrum (s. Eq. [TJ which reflects the flux 
emitted by a region which was not spatially resolved by the 
interferometer. An increasing effective baseline length of the 
interferometer results in a higher resolution. It has to be pointed 
out, that the single-dish spectra as well as the correlated spectra 
contain spectral contributions of the silicate emission from the 
whole disk, but the contributions from the hotter and brighter re- 
gions are increasing with increasing effective baseline length. In 
this context a homogeneous, axial-symmetric disk is assumed. 

According to the method used in Schegerer et al. (,200 6 ), we 
find a decreasing contribution of not-evolved, i.e. amorphous, 
0.1 ijm - small dust grains and an increasing crystallinity with 
increasing baseline length, i.e. decreasing distance to the cen- 
tral TTauri star (see Fig. |7] Table |5]). For comparison we add 
the single-dish, i.e. non-correlated spectrum Fy in Fig. [7] This 
single-dish spectrum confirms the derived tendencies. 

Figure |8] shows the crystallinity Cciys, which is plotted ver- 
sus the spatial resolution of our MIDI observations. The crys- 
talized material is concentrated mainly in the inner parts of the 
disk (point C for highest resolution), decreases strongly with de- 
creasing resolution (point B for the intermediate resolution) and 
converge to a lower limit for the single-dish observation (point 
A) which corresponds approximately to the abundance of crys- 
talline dust in interstellar matter (Gail |200"3T ). The relative mass 
contribution of small dust grains decreases from the outer to the 
inner disk regions. Considering the used spectroscopic slit we 
assume a resolution of 0.52 " for the single-dish observations. 
A corresponding result was previously found by van Boekel et 
al. (120041 ) for several HAeBe stars indicating more evolved sil- 
icate dust towards inner disk regions. Our result shows that the 
formation of crystalized dust grains is also favored in the in- 
nermost disk region of TTauri stars. A study of the absolute 
disk position of the crystalized dust, as for any other material, 
is out of the scope of this paper and should be presented in a 
future publication. Nonetheless such a forthcoming study is fa- 
vored by the fact that the position angles of our observations 
with MIDI are almost identical (Table |3]), i. e. Cci^s depends 
only on the radial coordinate r in the disk. Such a study would 
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Fig. 7. Single-dish (A) and correlated spectra (B) and (C) of the 
TTauri star RYTau (black, solid lines). The correlated spec- 
tra (B) and (C) corresponds to a baseline length of 49 m and 
78 m, respectively. The silicate feature are modeled by a lin- 
ear combination of mass absorption coefficients k, of different 
amorphous and crystalline silicates (grey lines). In order to ex- 
clude remnants of the data reduction we cut the wavelength in- 
terval ~9.0 fim to ~9.7 fim influenced by the telluric ozone band. 
Dashed lines and dot-dashed lines represent the contribution of 
0.1 jjm- and 1.5 /im-sized, amorphous grains, respectively. The 
dot-dot-dashed curves stand for the crystalline contribution. The 
underlying, solid curves represent the continuum (here: black- 
body with temperature T). 



Table 5. Results of our x^'^^ presented in Fig. Q The used 
method is described in Schegerer et al. ('2006) in detail. The un- 
derlying continuum is estimated by a single blackbody function 
with the temperature T. "RMC" stands for relative mass contri- 
bution, "am.-sma." for amorphous, 0.1 /im-small, "am.-la." for 
amorphous, 1.5 yum-large and "crys." for crystalline silicate dust 
grains. The crystalline component includes small and large sili- 
cate species: forsterite, enstatite, and quartz. See text for further 
discussion. 





baseline length 


resolution 


silicate 


RMC 


A 


single-dish 


70 AU 
(520 mas) 


am.-sma. 
am.-la. 
crys. 
T 


(15 + 3)% 
(82 ± 3)% 
(3 + 1)% 
(450 + 3) K 


B 


49 m 


2.8 AU 
(21 mas) 


am.-sma 
am.-la. 
crys. 
T 


(7 ± 3)% 
(82 ± 5)% 
(11 ±2)% 
(557 ± 5) K 


C 


78 m 


1.8 AU 
(13 mas) 


am.-sma 
am.-la. 
crys. 
T 


(1 ± 1)% 
(80 ± 4)% 
(19 ±2)% 
(722 + 9) K 
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Fig. 8. RMC of crystalized (dotted line, squares) and 0, 1 ;um- 
small, amorphous (dashed line, triangles) silicate grains plotted 
versus the reached spatial resolution of our observations. 



be an essential requirement to study the degree of radial mix- 
ing of material in circumstellar disks described by the ratio of 
the viscous inward and diffusive outwards stream (Wehrstedt 
& GailEOOll Bouw man et al. l2003l Gail|2004, Pavlyuchenko 
& Dullemond I2007i i. Moreover, such a study could give a hint 
whether circumstellar dust grains are only crystalized by thermal 
heating in the inner disk regions or whether electric discharges in 
outer disk regions can also crystalize dust grains. Both, a strong 
radial diffusion and electric discharges in outer disk regions, re- 
spectively, would imply a shallow decrease of the crystallinity 
with radius. Although Fig.|8]qualitatively shows that crystalized 
dust grains are located mainly in the inner disk region, a further 
study is required to determine the absolute disk position of the 
crystallised dust, along with more interferometric observations 
with different spatial resolutions. 

Finally, we add a few comments on the dust set used in our 
modeling approach. Instead of using the canonical dust set of 
"astronomical silicate" and carbon, we could also use the spe- 
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cific silicate dust composition found by the simple linear fit- 
ting routine presented above. Preibisch et al. (|20061 l have already 
used such an approach for modeling the HAeBe star HR5999. 
However, this proceeding is problematic in light of following ar- 
guments: 

i. The best-fit model parameters are not sensitive to very small 
grains (< 0.1 /vm). But such tiny grains are essential for 
generating the observed NIR emission as Weingartner & 
Draine (|2000l i pointed out. As mentioned in Sect. 14.21 a dust 
set where large dust particles a > I fim are implemented at 
the expense of tiny dust grains, is very ineffective in gener- 
ating NIR and MIR flux. Only additional modifications of 
the inner disk structure can compensate this lack of emis- 
sion. However, in disk systems where the NIR flux has al- 
ready vanished, the substition of small grains by large grains 
could be physically justified by advanced dust coagulation 
(Weinberger et al. 1999; Kornet et al. l200TT l. 

ii. The best-fit model parameters are not sensitive to carbon. 
The emission profile of carbon is strictly monotonic in the 
10 yum-wavelength range (e.g.. Wolf & Hillenbrand I2003I I 
and thus contributes to the underlying continuum, only. If 
carbon is considered in Eq. |7] too, each potential contribu- 
tion Ccarbon of carbon would not be independent of the con- 
tribution of the single blackbody function used to reproduce 
the underlying continuum. 

iii. Apart from the analysis of the contributions of each silicate 
component, another important result of our 10 //m-analysis 
is the derived tendency of increasing crystallinity towards 
inner disk regions. However, this result does not enclose the 
determination of the absolute disk position of the crystallised 
dust and it would be another requirement to use the results 
of our;i'^-fitting routine in our modeling approach. 

iv. It has not been clear, so far, if the determined contribution 
of each dust component corresponds to its mass fraction. As 
shown in Schegerer et al. (I2006I I the porosity of dust grains, 
which is not considered in our ;if^-fit, affects the shape of 
the silicate feature similar to the size effect. Therefore, the 
mass contributions of large, compact dust grains could be 
overestimated while the increasing crystallinity of the dust 
in the inner regions is still a safe conclusion (Sect.|6]l. 

Only further studies can clarify if and how the results of the pre- 
sented -fitting routine can be implemented in the modeling 
approach. 



7. Discussion 

The main aim of this paper was to model the structure of the cir- 
cumstellar environment of the T Tauri star RY Tau. In this con- 
text we presented different modeling approaches for the circum- 
stellar dust distribution (passive disk; active disk; active disk 
+ envelope) and pointed to potential supplements such as the 
puffed-up inner rim wall or the truncated outer disk model. An 
important aim of our approach was to keep the number of model 
parameters as small as possible. With respect to our modeling 
results additional parameters and modifications were only im- 
plemented if significant improvements could be obtained after- 
wards. In order to decrease the mm slope of the resulting SED, 
the disk with the MRN dust composition of astronomical silicate 
and carbon were replaced by a two-layer dust model where the 
disk interior also contains evolved, i.e. larger dust grains. 



7. 1 . The merits of the active disk model with and without an 
envelope. 

Accretion or dusty envelopes produce the additional infrared flux 
from A X 2 fim up to /i ^ 8 fim of the SED that is missed in the 
model of a "naked", passively heated disk. 

The extra infrared radiation from the active disk model is 
generated in the innermost disk region inwards the point of dust 
sublimation (AppendixlAb. The radiation which is caused by ac- 
cretion additionally heats the innermost disk regions close to the 
inner disk radius. In comparison to a "naked", passively heated 
disk, the implementation of accretion effects in the model re- 
sults in a higher spatial concentration of infrared flux in the in- 
ner region that is not spatially resolved by our MIDI measure- 
ment. The computed MIR visibility increases, therefore (com- 
pare middle and lower panel in Fig. |6ll. The accretion rate 
of the model presented here (active disk without envelope), 
9.1 X lO^^Moyr"', is smaller than the value found by Akeson 
et al. ( 120051 2.5 x lO^^Moyr"^) in their model study but corre- 
sponds to that of Calvet et al. (l2D04l 6.4-9. 1 x 10 ^^Moyr '). The 
latter result is based on a multi-wavelength study in the optical- 
UV range considering different emission-line profiles. An accre- 
tion rate of 2.5 x lO^^Moyr"' decreases the visibilities in the 
MIR range because of a stronger irradiation of outer disk re- 
gions. Based on a flux ratio measurement between the contin- 
uum excess and the intrinsic photospheric flux at a wavelength 
of 5700 A, Hartigan et al. ( 1995) found a much lower accretion 
rate of M ^ 2.5 x lO"** Moyr"' in RYTau. Vink et al. ( |2003] ) 
determined an accretion rate of M a: 7.5 x 10"*^ Moyr~'. The 
derived accretion rate from the model of Akeson et al. 120051 is 
up to a factor of 10 larger than the values that were measured 
in the UV range. Such a discrepancy between standard accretion 
disk models and the measurements has been already discussed 
by Muzerolle et al. (120041 ). To obtain a consistency in the mea- 
surements, they introduced an artificially puffed-up inner rim in 
their modeling approach, accounting for the large NIR excesses 
of classical T Tauri stars, but without requiring excessive accre- 
tion rates. Instead of implementing such a puffed-up inner rim 
we assumed an envelope in our active disk -i- envelope where 
we assumed a smaller accretion rate (factor ~4) than in the pure 
active disk model. 

NIR flux in the model with the envelope has its origin close 
to but outwards of the sublimation point. This extra flux results 
from the stellar heating of the dust in the envelope up to the 
sublimation temperature and in an increase of the MIR visibility 
corresponding to the pure active model. 

In contrast to pure accretion, a dusty envelope dims the cen- 
tral star and prevents the outer disk regions from being heated 
too strongly by direct stellar irradiation. Therefore, the MIR ree- 
mission from these outer regions is decreased and the spatial 
concentration of the infrared radiation in the inner regions is in- 
creased. Such an effect caused by the envelope results in an in- 
crease of the MIR visibility, in particular for measurements with 
the smaller projected baseline of 48 m. We have to mention that 
comparable effects could also be achieved by the truncation of 
the outer disk in our active disk model or by a strongly puffed-up 
inner rim. 

In contrast to our finding Akeson et al. ( 120051 ) did not find 
any hints for an additional envelope in their modeling approach 
for RYTau based on NIR, interferometric observations. They 
suggest NIR unveiled CO absorption lines which RY Tau is ex- 
hibiting. In fact, such unveiled CO lines could be evidence of the 
absence of a substantial envelope (Najita et al. 2003; Calvet et 
al. il997i l. However, the almost unveiled CO absorption lines do 
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not necessarily exclude an envelope when considering the fol- 
lowing argument: 

The results of two-dimensional collapse calculations of the 
infalling matter in an envelope (Yorke et al. 1993) suggests a 
more plane-parallel than spherical envelope geometry. In this 
model, the envelope has already collapsed at the inner disk edge. 
With MIDI we resolved inner disk regions with a distance of 
several AU from the star (see Fig. |2]) where remnants of the 
envelope could still exist. However, at these distances the gas 
in the potential envelope is already too cold to provide a sub- 
stantial veiling of CO lines. We point to a study of Bastien & 
Landstreet ( 1979) where it was suggested that most of the polar- 
ization found towards RY Tau arises from a circumstellar (dusty) 
envelope which actually lies outside of the high-temperature, 
gas-emitting region. Therefore, a geometry where the envelope 
has already disappeared at the inner edge but not at adjacent re- 
gions could explain why the observed NIR CO absorption lines 
are unveiled and why Akeson et al. (.2005) failed to model pho- 
tometric and NIR visibility data considering an additional en- 
velope to their pure active disk model. In fact, with respect to 
Fig.|2]the NIR emission mainly originates close to the sublima- 
tion point. 

Another potential origin of such an (not necessarily spheri- 
cal) envelope could be magnetically driven disk winds contain- 
ing gas and, additionally, small quantities of small dust parti- 
cles: material at the disk surface could even follow magnetical 
field lines for radii r which are much larger than the magneti- 
cal truncation radius r » 7?bnd (Appendix lAl) and larger than the 
inner radius r > Ri„ (Blandford & Payne 1982 ). Former studies 
assumed a correlation between the disk accretion rate M onto 
the star and the outflow mass-loss rate Mof with Mof/M ~ 0.1 
(Richer et al. 120001 ). These dust particles that follow the outflow- 
ing wind from the disk form the optically thin dusty envelope 
assumed in our modeling approach. A possible consequence of 
this procedure would be the (acceleration of the) formation of 
an inner gap in the innermost disk region as observed in older 
TTauri objects such as TWHya (Calvet et al. 120021 ). Finally, 
we note that Fendt & Camenzind (1996) have studied stationary, 
axisymmetric wind flows driven by a rapidly rotating magneto- 
sphere. They found that the (gas) particle density in the outflow 
decreases with r^^-^. This result is independent from the stellar 
parameters (private communication with Ch. Fendt). 

Another origin of the dust particles in the circumstellar enve- 
lope could be hard UV irradiation from the inner accretion zone 
or/and the star which increases the gas temperature in the up- 
per layers of the disk up to ~10'* K allowing the hot, gaseous 
material to escape from the gravitationally bound system of the 
star Small dust particles accompany the gas outflow. This effect, 
called photoevaporation, effectively starts at a critical radius 
where the sound speed is in the range of the escape speed, i. e. 
a few AU in the case of RYTau (see Dullemond et al. '2006' for 
a review). However, only further theoretical studies can clear- 
ify if the mass loss rates of dusty material in these disk out- 
flows is high enough to cause a sufficiently strong effect visible 
with MIDI. For completeness we note that gas pressure depen- 
dent photophoretic forces of light can also induce the ejection of 
dust from the optical thin surface layers of the disk as studied by 
Wurm (120051 ) and Wurm & Krauss (2006). 

YSOs of Class I (Adams et al. 1987) that reveal a circumstel- 
lar envelope + disk structure typically also show signs for accre- 
tion given that accretion is stronger for younger objects (e.g., 
Hartmann et al. |2005| ). Accretion is certainly present in RY Tau, 
too, according to the mentioned, numerous studies, such as the 
analysis of the existing Bry and Halpha lines and UV excess 



radiation. However, we have also noted that the implementa- 
tion of accretion can be totally ignored in an envelope + pas- 
sive disk model to reproduce SED and MIR visibilities. Both, 
accretion and envelope, increase the NIR- and MIR-flux. Only 
complementary observations in the UV range where the accre- 
tion rate can be independently measured, will provide additional 
constraints to disentangle the different model approaches and al- 
low us to consider both accretion and the envelope in one model. 

7.2. A potential Stellar companion? 

As mentioned in Sect. |2] Bertout et al. ( 119991 ) found indi- 
rect hints for a stellar companion analysing HIPR\RCOS data. 
Assuming the regular motion of the photocenter of RY Tau they 
derived a projected minimum distance of 3.27 AU and a position 
angle of 304° + 34° for the potential secondary. The method they 
used is described in Wielen (119961 ). 

A potential detection of a companion by interferometric ob- 
servations depends on the separation Osep, the position angle of 
the companion with respect to the position angle of the interfer- 
ometric baseline and the brightness ratio Oiat of secondary and 
primary. The visibility of a binary system can be expressed by 
the approximate formula: 



J 1 + fl^^j + 2arat COS(27rasepf ) 

V(B, A) = ao(B, A)^ . (8) 

1 + fli-at 

The parameter ao(B, A) represents the visibility that results from 
the circumstellar and -binary material assumed to be equally 
distributed around the components. A positive detection could 
then be recognized by sinusoidal variation of the visibility 
(Ratzkal 2005b . Furthermore, all interferometric observations are 
snapshots of the system in contrast to the HIPPARCOS observa- 
tions repeated on long-term. A disadvantageous configuration of 
the secondary could actually prevent detection with the interfer- 
ometer. With the non-dection of a clear binary signal in our MIDI 
data, we can thus neither verify, nor disprove the existence of a 
companion. Assuming flo(B, A) = I for two point sources, Fig.|9] 
shows the relation between the interferometric binary signal and 
the projected baseline and the separation of the objects, respec- 
tively. With increasing baseline length or increasing separation 
the frequency of the sinusoidal variation also increases. The am- 
plitude of the variation decreases with the brightness ratio of the 
components. 

Another aspect of this discussion is the possibility that the 
regular motion of the photocenter observed with HIPPARCOS 
in the visible range is not caused by a secondary, but by a bright- 
ness irregularity in the circumstellar environment, i.e. in the disk 
and/or envelope. In fact, considering scattered-light images of 
Class I objects (Padgett et al. 1 999) , circumstellar envelopes, 
which have their brightness maximum in the NIR and adjacent 
wavelength ranges (Wolf et al. 2003), reveal such brightness 
irregularities. However, it is not clear if such a brightness ir- 
regularity in the circumstellar environment of RY Tau is strong 
enough to cause the observed regular motion of the photocenter. 

7.3. Comparison with HAeBe stars 

HAeBe stars were formerly classified in Group I and Group 
II sou rces (M eeus et al. '2001; Sect. 15^ . Recendy, Leinert 
et al. ( |2004| ) could explain this phenomenological classifica- 
tion after modeling the interferometric measurements of several 
HAeBe objects obtained with MIDI. In fact, they noticed that 
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Fig. 9. Theoretical prediction for the result obtained for the ob- 
servation of a stellar binary sytem with MIDI. Here, we assume 
two point sources, i. e. ao(B,A) = 1 in Eq.[8] The (normalized) 
visibility is shown in relation with the separation of both stars 
(solid lines: Osep - 4 AE; dashed lines: Osep = 10 AE) and the 
projected baseline length (left: B = 49 m; right: B - 79 m). The 
brightness ratio of the components are 1 : 3. In this example the 
baseline of the interferometer is also parallel to the connection 
line of the components. 



the half light radius of their disk models in the MIR wavelength 
rang^ linearly correlates with the IRAS color between 12 yum 
and 25 yum, i. e. -2.5 log(Fv(12 jum)/Fy(25 fim)). Circumstellar 
disks around Group II sources are less flared. The geometrical 
effect results in the outer disk regions around Group II sources 
being less strongly heated than the corresponding disk regions of 
Group I sources. Therefore, the size of the half light radius and, 
simultaneously, the MIR color is smaller for Group II sources. 

Because of its stellar luminosity (L* = 11.5 Lq) and stellar 
mass (Mi, = 1.69 Mq) the TTauri star RYTau could be con- 
sidered as a transition object between TTauri and HAeBe stars. 
As the FIR-flux declines, RYTau can be classified as a Group 
II source (Meeus et al. |200TT l. The IRAS color of RYTau in 
the MIR wavelength range is -2.5 log{Fy(l2 fj.m)/Fy{25 jt/m)) - 
0,43 + 0, 11 (IRAS Catalogs 1985). The half light radius th, i.e. 
the FWHM of the intensity distribution in our best active disk 
model is determined here by a Gaussian fit. The intensity dis- 
tribution is computed for a wavelength of A - 12.5 fim as it 
is not affected by the silicate emission at this wavelength. We 
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* The half light radius encircles the disk region emitting half of the 
totally released energy in the disk. 



Fig. 10. Correlation between the half light radius ru at a wave- 
length of 12.5 jjm and the IRAS-color at /i = 12 fim and 
A - 25 /vm. Plotted are the corresponding values for HAeBe 
objects (Leinert et al. 1 2004] ) and for RY Tau. The error bars cor- 
respond to errors of the IRAS fluxes and of the interferometric 
measurements. The value ru of RYTau is an upper limit, as the 
active disk model predicts a too small visibility in comparison to 
the measurement. 



obtain FWHM - 1 .75 AU. This value is an upper limit of the 
half light radius as the active disk model predicts a too small 
visibility in comparison to the measurement, i.e. the intensity 
distribution of the model disk decreases too slowly for increas- 
ing radii in contrast to the real intensity distribution (measure- 
ment). However, with respect to the study of Leinert et al. ( 120041 
Fig. 5) this result actually confirms the linear correlation for- 
merly found for HAeBe stars, solely. A further correspondence 
between T Tauri and HAeBe objects has been already shown in 
Schegerer et al. ( I2006I I considering correlations between stellar 
properties and the silicate composition of the circumstellar disk 
of TTauri and HAeBe objects. 

7.4. Complex interplay of disk parameters 

The visibility, which we try to model, is a complex function of 
many different disk and dust parameters. In such complex disk 
models the effect of any modifications of geometrical disk pa- 
rameters cannot always be predicted. However, the tendencies 
of modifications of the main model parameters are summarized 
in the following: 

i. An increase of the stellar effective temperature and stel- 
lar luminosity L*, in particular, mainly causes an increase 
of the visible and infrared flux. If the flux of the unresolved 
star increases more strongly than the heating of the dust and 
the MIR flux from the resolved circumstellar surroundings, 
the MIR-visibility increases. The accretion luminosity af- 
fects SED and MIR visibility, correspondingly. An increase 
of the stellar mass corresponds to an increase of the ac- 
cretion luminosity in the active disk model ( Eq. lA.2b . 

ii. The optical depth of the disk increases with the disk mass 
assuming a constant outer radius. Compacter disk structure 
decreases the averaged temperature in the inner disk region 
where less MIR flux is emitted. However, the MIR flux from 
outer regions simultaneously increases and the MIR visibil- 
ity therefore decreases. 
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iii. A decrease of the inner disk radius Rm results in more NIR 
and MIR flux and slightly increases the MIR visibility. 

iv. As shown in Sect. 15.11 the outer disk radius determines the 
FIR flux in our modeling approach. The eff'ect of an increase 
of the outer radius corresponds to the effect of a decrease of 
the disk mass. 

V. Similar to the outer disk radius, the exponent p has distinct 
effects on the radial density distribution in the disk. Surface 
density distributions 2 with larger gradients p result in an 
increase of the density in the inner regions of the disk. The 
MIR visibility and the infrared excess from the disk in the 
NIR up to the MIR wavelength range increases, simultane- 
ously. While the FIR excess can decrease strongly at values 
greater than p - 2.0, smaller exponents p < 0.8 only negli- 
gibly affect the SED and MIR visibility. 

The spatial information, which we got from the MIDI observa- 
tions, strongly reduced the number of disk models which can 
reproduce the SED. Certainly, the outer disk radius and the disk 
mass are less constrained parameters in our models. The unique- 
ness of our approach cannot be proven. However, in our ap- 
proach we consistently verified whether the modeling results can 
be improved by varying the modeling parameters. The step sizes 
used for a final variation of these modeUng parameters are: 

A/?out = 10 AU, ARin = 0.05 AU, 

AL = 0.1 Lq, Ac2 = 0.5, and Ap = 0.1. 

For the parameter ci and Mjisk we used: 

Aci = O.lci, and AMdisk - 0.5Mdisk, 

respectively. These step sizes can be considered as an error of 
each resulting parameter (Table|4|i. As mentioned above, the stel- 
lar mass Mstar, the effective stellar temperature T*, the accretion 
rate M, the boundary temperature Thnd and radius /?bnd as well 
as the visual extinction Ay are constant parameters considering 
previous measurements. 

8. Summary 

We present interferometric observations of the classical T Tauri 
star RY Tau in the 10 fim range, which show the source well re- 
solved, together with the total spectrum of the source. We modi- 
fied the MC3D code (Wolf et al. 1999) to obtain a self-consistent 
model (in the temperature and density distribution) of the cir- 
cumstellar disk around RY Tau, including accretion, with the fol- 
lowing results: 

i. Both, the SED and the 10 /urn interferometric measurements 
can be fitted by a self-consistent, physically reasonable ra- 
diative transfer model. The accretion and an envelope have 
the effect of filling in the missing emission in the range from 
3 fim up to 8 /im. 

ii. A model of a circumstellar active disk with an optically thin 
envelope and a low accretion rate represents the data and the 
model without the envelope and a much larger accretion rate 
(Akeson et al. 2005). The modeling approach with an enve- 
lope, which could have its origin in an outflowing gaseous 
and dusty disk wind, should be considered as a future, viable 
alternative for the puffed-up inner rim, which results from 
a vertical expansion of the disk edge according to Natta et 
al. (1200 11 1 and Dullemond et al. (2001). The latter approach 
could not be confirmed with our particular model of RY Tau 
considering a self-consistently determined density distribu- 
tion considering hydrostatical equilibrium. 



iii. The interferometric observation reduces the number of mod- 
els that solely reproduce the SED of RY Tau. 

iv. The results of our modeling approaches are: The active 
disk with or without an (optically thin) envelope predicts 
a system with an upper limit of ^ 55° for the inclina- 
tion. While a stellar mass of M* - 1.69 Mq and an ef- 
fective temperature of 5560 K can be adopted from previ- 
ous measurements (Table [TJ, we found a stellar luminos- 
ity of Li, w 10. L0 which is smaller than the estimates of 
Akeson et al. (120051 12.8 Lq) but corresponds to the find- 
ings of Calvet et al. (120041 9.6 L©). The inner radius of our 
models Rin - 0.3 AU coincides with the result of Akeson 
et al. (2005, Rin = 0.27 AU) where interferometrical studies 
in the NIR preceded. The accretion rate of the model with- 
out an envelope is 9.1 x lO^^Moyr"' and corresponds to the 
upper limit found by Calvet et al. (|2004| l. The use of an ad- 
ditional envelope allows us to use a lower accretion rate of 
2.5 X lO^'^Moyr"' which is the same as the result published 
by Hartigan et al. ( I1995I I. 

V. Comparison of interferometric and single-dish observations 
shows for the first time dust evolution in a T Tauri star with a 
reduced fraction of small amorphous and an increased frac- 
tion of crystalline particles closer to the star. This is similar 
to the finding for HAeBe stars (van Boekel et al. 120041 ) . 

vi. Our data neither support nor contradict the existence of a 
close (~3 AU) companion for which some astrometric evi- 
dence exists (Bertout et al. |19991 ). 

vii. The complex interplay of different disk parameters is con- 
sistently apparent from the use of the Monte Carlo radiative 
transfer code. 
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Appendix A: The accretion model 

Accretion models have been studied in many publications and 
the parameters that we used for our active disk model are well- 
known properties of such models. However, only a few authors 
have implemented accretion effects in radiative transfer models 
(e.g., Akeson et al. 2005 ). It is generally assumed that the pas- 
sive heating of the disk is the dominating source of infrared and 
mm irradiation. This is certainly true for the large scales of the 
disk but not for its innermost region (< 2AU; e.g., D'Alessio et 
al.[1998 ) which can be resolved with long-baseline interferome- 
ters. We briefly summarize here the parameters of the accretion 
disk model implemented in our approach. 

Former accretion disk models are generally based on the 
assumption of a geometrically thin, steady disk established by 
Lynden-Bell & Pringel ( 119741 ) and Pringle ( il9811 l. In this canon- 
ical accretion model it is assumed that viscous stresses within the 
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disk transport angular momentum to its outer regions. As a con- 
sequence of this, most of the disk material moves inward onto 
the protostar, while some disk matter moves outward, absorbing 
all the angular momentum. Assuming a geometrically thin and 
steady disk the conservation of transversal and angular momen- 
tum results in the dissipation rate D per unit area and time as a 
function of the radial distance from the star r. 



D(r) 



3GM^M 
Anr^ 



1 



(^) 



(A.l) 



The quantities G, M*, M, r and 7?* are in this order: gravity con- 
stant, stellar mass, accretion rate of the infalling material, radial 
distance and stellar radius. The total released energy between an 
inner boundary radius of the disk /?bnd and infinity is 



■^'disk = I 



D{r) 2nrdr — 



1 GM^M 

2 Rbnd 



(A.2) 



Lynden-Bell & Pringle ( 119741 1 assumed that the disk extends 
down to the stellar radius /?bnd = R*, or to the corotation ra- 
dius where the Keplerian orbital period equals the stellar rotation 
period. Recent observations (e.g., Muzerolle et al. 2004) con- 
firmed a magnetically mediated accretion processes in the inner- 
most disk. The accretion disk is truncated by the stellar magnetic 
field at a radius /?bnd of several stellar radii and material is chan- 
neled along magnetic field lines onto the star with nearly free- 
fall velocity (e.g., Uchida & Shibata 1984; Bertout et al. 1988)1. 
The magnetic flux density B of a T Tauri star is in the range of 
~1 - 3 kG (e.g., Johns-Krull et al. I2003I I. As the boundary ra- 
dius /?bnd is assumed to be smaller than the sublimation radius 
of dust {~Rin for RYTau), the accreting material in this inner 
region is gaseous. Close to the stellar surface the accreting, free- 
falling material is abruptly stopped in one or a series of shock 
layers. In this inner accretion zone the total released energy is 
^GM^M/R^ with ^ = 1 - R-^/{2Rbnd)- Therefore, adjacent lay- 
ers of the shock region like the stellar photosphere are heated to 
temperatures of Tbnd - 5700 K up to 8800 K, while tempera- 
tures of 9000 up to 20, 000 K are achieved in upper shock areas 
(Calvet & Gullbring |Tg98] l. Muzerolle et al. ( |20D4] l found in their 
studies that the energy contribution from the upper shock layers 
can be neglected for accretion rates M ^ 10"^ Moyr"' while the 
stellar photosphere emits most of the accretion energy. Thus, the 
stellar photosphere is assumed to emit the intrinsic stellar flux 
and the accretion energy (GMi,M /R-,,. 

The crux of the accretion theory is the surface coverage fac- 
tor / which mimics an annulus on the stellar photosphere where 
the additional accretion emission occurs (D'Alessio et al. |1998l 
Bertout et al. 1988 ; Lynden-Bell & Pringle [19741 1. If a blackbody 
emitter is assumed, one gets: 



f4nRlaTl 



(A.3) 



The Stefan-Boltzmann constant is represented by cr. Calvet & 
Gullbring (199E) determined / « 0.1% - 1.0% for HAeBe stai's 
and / ^ 10% for T Tauri stars. Stronger veiling of absorption 
lines indicates the larger coverage factors / of T Tauri stars. 

In the region between the sublimation and magnetic bound- 
ary radius the gradually released accretion rate in our model (see 
Eq. I A. Il l is determined by the accretion theory of Lynden-Bell 
& Pringles (1974) . Because of the steep decrease of the dissi- 
pation rate D with increasing distance r (see Eq. lA.lb accretion 
at radial distances r > Rin is neglected in our model. The re- 
maining gravitational energy of the accreting material between 
r - Rbnd and r - Ri, is emitted at the stellar surface assum- 
ing a blackbody emitter with the temperature T^nd- Note that the 



effects of the boundary temperatur Tbnd and the magnetic trun- 
cation radius /?bnd on the SED and MIR visibilities are marginal 
(see D'Alessio et al. |1998l l. A sketch of our active disk model is 
shown in Fig. lA.ll 
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Table 4. Parameter set for the passive disk model, the active model ("Sect. |5TT] ) and the active model with an envelope ("Sect. l53T l. 
Results correspond to Fig.|6] 



parameter 


passive disk 


active disk 


active disk + envelope 




model 


model 


model 


Stellar mass Mj, 


1.69 Mo 


1.69 Mo 


1.69 Mo 


Stellar temperature Ty, 


5560 K 


5560 K 


5560 K 


Stellar luminosity Ly, 


10.0 Lo 


10.0 Lo 


10.0 Lo 


Disk mass Mjjsi; 


1 X IQ-- Mo 


4 X lO-'' Mo 


4 X 10-3 Mo 


Outer disk radius R„ut 


270 AU 


270 AU 


270 AU 


Inner disk radius R^y, 


0.3 AU 


0.3 AU 


0.3 AU 


Exponent p (see Eq.|2ll 


1.3 


1.3 


1.4 


Inclination & 


g 70° 


S 70° 


S65° 


ci (see Eq.|5]( 






5.0 X 10-5 


C2 (see Eq.|5ll 






1.0 


Accretion rate M 




9.1 X 10-**Moyr-' 


2.5 X 10-'*Moyr"' 


Boundary temperature Tbnd 




8000 K 


8000 K 


Truncation radius R^nd 




5R* 


5R* 




Fig. 6. Top row: Models of the spectral density distribution of RYTau for an inclination angle of - 25°. Here, we assume a 
"naked", passively heated disk model without considering accretion effects and an envelope. Real photometric data with error bars 
are included in both models (see Table|2]i. The dashed line represents the intrinsic, stellar flux. The flux F,i is given in units Wm"^. 
According to the theorem of van Cittert-Zernicke the modeled visibilities in all diagrams were calculated from the corresponding 
model image for the wavelengths of 8.5, 9.5, 10.6, 11. 5 and 12.5 /urn. Triangles and squares represent the upper and lower limit of 
the visibilities V(A) for different position angles but the same inclination of the model. The measured data are included with error 
bars. 

Middle row: Model of the spectral density distribution of RY Tau for an inclination of 25° considering accretion effects in addition 
to a passively heated disk. 

Bottom row: SED that is obtained only from an active disk model considering an envelope, additionally. 
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gas & disk gas & dust 




interior disk layers 

Fig. A.l. Simplified sketch of our active disk model. The boundary temperature Tbnd, the boundary radius 7?bnd and the inner disk 
radius Ri„ are shown. We use a two-layer disk model (Sect. l4!2] ): a disk atmosphere where the optical depth tn is smaller than unity 
in N band and a optical thick, disk interior. 



